Neutralisation sensitivity of the SARS-CoV-2 BA.2.87.1 variant (preprint)

Against the backdrop of the rapid global takeover and dominance of BA.1/BA.2 and subsequently BA.2.86 lineages, the emergence of a highly divergent SARS-CoV-2 variant warrants characterization and close monitoring. Recently, another such BA.2 descendent, designated BA.2.87.1, was detected in South Africa. Here, we show using spike-pseudotyped viruses that BA.2.87.1 is less resistant to neutralisation by prevailing antibody responses in Sweden than other currently circulating variants such as JN.1. Further we show that a monovalent XBB.1.5-adapted booster enhanced neutralising antibody titers to BA.2.87.1 by almost 4-fold. While BA.2.87.1 may not outcompete other currently-circulating lineages, the repeated emergence and transmission of highly diverged variants suggests that another large antigenic shift, similar to the replacement by Omicron, may be likely in the future.

Molecular Epidemiology and Diversity of SARS-CoV-2 in Ethiopia, 2020-2022 (preprint)

Ethiopia is the second most populous country in Africa and the sixth most affected by COVID-19 on the continent. Despite having experienced five infection waves, >499 000 cases, and ~7 500 COVID-19-related deaths as of January 2023, there is still no detailed genomic epidemiological report on the introduction and spread of SARS-CoV-2 in Ethiopia. In this study, we reconstructed and elucidated the COVID-19 epidemic dynamics. Specifically, we investigated the introduction, local transmission, ongoing evolution, and spread of SARS-CoV-2 during the first four infection waves using 353 high-quality near-whole genomes sampled in Ethiopia. Our results show that whereas viral introductions seeded the first wave, subsequent waves were seeded by local transmission. The B.1.480 lineage emerged in the first wave and notably remained in circulation even after the emergence of the Alpha variant. The B.1.480 was out-competed by the Delta variant. Notably, Ethiopia lack of local sequencing capacity was further limited by sporadic, uneven, and insufficient sampling that limited the incorporation of genomic epidemiology in the epidemic public health response in Ethiopia. These results highlight Ethiopia role in SARS-CoV-2 dissemination and the urgent need for balanced, near-real-time genomic sequencing.

Global Expansion of SARS-CoV-2 Variants of Concern: Dispersal Patterns and Influence of Air Travel (preprint)

In many regions of the world, the Alpha, Beta and Gamma SARS-CoV-2 Variants of Concern (VOCs) co-circulated during 2020-21 and fueled waves of infections. During 2021, these variants were almost completely displaced by the Delta variant, causing a third wave of infections worldwide. This phenomenon of global viral lineage displacement was observed again in late 2021, when the Omicron variant disseminated globally. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of SARS-CoV-2 VOCs worldwide. We find that the source-sink dynamics of SARS-CoV-2 varied substantially by VOC, and identify countries that acted as global hubs of variant dissemination, while other countries became regional contributors to the export of specific variants. We demonstrate a declining role of presumed origin countries of VOCs to their global dispersal: we estimate that India contributed <15% of all global exports of Delta to other countries and South Africa <1-2% of all global Omicron exports globally. We further estimate that >80 countries had received introductions of Omicron BA.1 100 days after its inferred date of emergence, compared to just over 25 countries for the Alpha variant. This increased speed of global dissemination was associated with a rebound in air travel volume prior to Omicron emergence in addition to the higher transmissibility of Omicron relative to Alpha. Our study highlights the importance of global and regional hubs in VOC dispersal, and the speed at which highly transmissible variants disseminate through these hubs, even before their detection and characterization through genomic surveillance.

Variable loss of antibody potency against SARS-CoV-2 B.1.1.529 (Omicron) (preprint)

The recently-emerged SARS-CoV-2 B.1.1.529 variant (Omicron) is spreading rapidly in many countries, with a spike that is highly diverged from the pandemic founder, raising fears that it may evade neutralizing antibody responses. We cloned the Omicron spike from a diagnostic sample which allowed us to rapidly establish an Omicron pseudotyped virus neutralization assay, sharing initial neutralization results only 13 days after the variant was first reported to the WHO, 8 days after receiving the sample. Here we show that Omicron is substantially resistant to neutralization by several monoclonal antibodies that form part of clinical cocktails. Further, we find neutralizing antibody responses in pooled reference sera sampled shortly after infection or vaccination are substantially less potent against Omicron, with neutralizing antibody titers reduced by up to 45 fold compared to those for the pandemic founder. Similarly, in a cohort of convalescent sera prior to vaccination, neutralization of Omicron was low to undetectable. However, in recent samples from two cohorts from Stockholm, Sweden, antibody responses capable of cross-neutralizing Omicron were prevalent. Sera from infected-then-vaccinated healthcare workers exhibited robust cross-neutralization of Omicron, with an average potency reduction of only 5-fold relative to the pandemic founder variant, and some donors showing no loss at all. A similar pattern was observed in randomly sampled recent blood donors, with an average 7-fold loss of potency. Both cohorts showed substantial between-donor heterogeneity in their ability to neutralize Omicron. Together, these data highlight the extensive but incomplete evasion of neutralizing antibody responses by the Omicron variant, and suggest that increasing the magnitude of neutralizing antibody responses by boosting with unmodified vaccines may suffice to raise titers to levels that are protective.

Conserved recombination patterns across coronavirus subgenera (preprint)

Recombination contributes to the genetic diversity found in coronaviruses and is known to be a prominent mechanism whereby they evolve. It is apparent, both from controlled experiments and in genome sequences sampled from nature, that patterns of recombination in coronaviruses are non-random and that this is likely attributable to a combination of sequence features that favour the occurrence of recombination breakpoints at specific genomic sites, and selection disfavouring the survival of recombinants within which favourable intra-genome interactions have been disrupted. Here we leverage available whole-genome sequence data for six coronavirus subgenera to identify specific patterns of recombination that are conserved between multiple subgenera and then identify the likely factors that underlie these conserved patterns. Specifically, we confirm the non-randomness of recombination breakpoints across all six tested coronavirus subgenera, locate conserved recombination hot- and cold-spots, and determine that the locations of transcriptional regulatory sequences are likely major determinants of conserved recombination breakpoint hot-spot locations. We find that while the locations of recombination breakpoints are not uniformly associated with degrees of nucleotide sequence conservation, they display significant tendencies in multiple coronavirus subgenera to occur in low guanine-cytosine content genome regions, in non-coding regions, at the edges of genes, and at sites within the Spike gene that are predicted to be minimally disruptive of Spike protein folding. While it is apparent that sequence features such as transcriptional regulatory sequences are likely major determinants of where the template-switching events that yield recombination breakpoints most commonly occur, it is evident that selection against misfolded recombinant proteins also strongly impacts observable recombination breakpoint distributions in coronavirus genomes sampled from nature.

Rapid replacement of the Beta variant by the Delta variant in South Africa (preprint)

The Beta variant of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in South Africa in late 2020 and rapidly became the dominant variant, causing over 95% of infections in the country during and after the second epidemic wave. Here we show rapid replacement of the Beta variant by the Delta variant, a highly transmissible variant of concern (VOC) that emerged in India and subsequently spread around the world. The Delta variant was imported to South Africa primarily from India, spread rapidly in large monophyletic clusters to all provinces, and became dominant within three months of introduction. This was associated with a resurgence in community transmission, leading to a third wave which was associated with a high number of deaths. We estimated a growth advantage for the Delta variant in South Africa of 0.089 (95% confidence interval [CI] 0.084-0.093) per day which corresponds to a transmission advantage of 46% (95% CI 44-48) compared to the Beta variant. These data provide additional support for the increased transmissibility of the Delta variant relative to other VOC and highlight how dynamic shifts in the distribution of variants contribute to the ongoing public health threat.

A Novel and Expanding SARS-CoV-2 Variant, B.1.1.318, dominates infections in Mauritius (preprint)

Mauritius, a small island in the Indian Ocean, has had a unique experience of the SARS-CoV-2 pandemic. In March 2020, Mauritius endured a small first wave and quickly implemented control measures which allowed elimination of local transmission of SARS-CoV-2. When borders to the island reopened, it was accompanied by mandatory quarantine and testing of incoming passengers to avoid reintroduction of the virus into the community. As variants of concern (VOCs) emerged elsewhere in the world, Mauritius began using genomic surveillance to keep track of quarantined cases of these variants. In March 2021, another local outbreak occurred, and sequencing was used to investigate this new wave of local infections. Here, we analyze 154 SARS-CoV-2 viral genomes from Mauritius, which represent 12% of all the infections seem in Mauritius, these were both from specimens of incoming passengers before March 2021 and those of cases during the second wave. Our findings indicate that despite the presence of known VOCs Beta (B.1.351) and Alpha (B.1.1.7) among quarantined passengers, the second wave of local SARS-CoV-2 infections in Mauritius was caused by a single introduction and dominant circulation of the B.1.1.318 virus. The B.1.1.318 variant is characterized by fourteen non-synonymous mutations in the S-gene, with five encoded amino acid substitutions (T95I, E484K, D614G, P681H, D796H) and one deletion (Y144del) in the Spike glycoprotein. This variant seems to be increasing in prevalence and it is now present in 34 countries. This study highlights that despite having stopped the introduction of more transmissible VOCs by travel quarantines, a single undetected introduction of a B.1.1.318 lineage virus was enough to initiate a large local outbreak in Mauritius and demonstrated the need for continuous genomic surveillance to fully inform public health decisions.

A novel variant of interest of SARS-CoV-2 with multiple spike mutations is identified from travel surveillance in Africa (preprint)

At the end of 2020, the Network for Genomic Surveillance in South Africa (NGS-SA) detected a SARS-CoV-2 variant of concern (VOC) in South Africa (501Y.V2 or PANGO lineage B.1.351)1. 501Y.V2 is associated with increased transmissibility and resistance to neutralizing antibodies elicited by natural infection and vaccination2,3. 501Y.V2 has since spread to over 50 countries around the world and has contributed to a significant resurgence of the epidemic in southern Africa. In order to rapidly characterize the spread of this and other emerging VOCs and variants of interest (VOIs), NGS-SA partnered with the Africa Centres for Disease Control and Prevention and the African Society of Laboratory Medicine through the Africa Pathogen Genomics Initiative to strengthen SARS-CoV-2 genomic surveillance across the region. Here, we report the first genomic surveillance results from Angola, which has had 21 500 reported cases and around 500 deaths from COVID-19 up to March 2021 (Supplemental Fig S1). On 15 January 2021, in response to the international spread of VOCs, the government instituted compulsory rapid antigen testing of all passengers arriving at the main international airport, in addition to the existing requirement to present a negative PCR test taken within 72 hours of travel. All individuals with a positive antigen test are isolated in a government facility for a minimum of 14 days and require two negative RT-PCR tests at least 48 hours apart for de-isolation, whilst all travelers with a negative test on arrival proceed to mandatory self-quarantine for 10 days followed by a repeat test. In March 2021, we received 118 nasopharyngeal swab samples collected between June 2020 and February 2021, a number of which were from incoming air travelers (Supplemental Fig S1). From these, we produced 73 high quality genomes (>80% coverage), 14 of which were known VOCs/VOIs (seven 501Y.V2/B.1.351, six B.1.1.7, one B.1.525), 44 of which were C.16 (a common lineage circulating in Portugal), and twelve of which were other lineages (Supplemental Fig S2). In addition, we detected a new VOI in three incoming travelers from Tanzania who were tested together at the airport in mid-February. The three genomes from these passengers were almost identical and presented highly divergent sequences within the A lineage (Figure 1A & 1B). The GISAID database contains nine other sequences reported to be sampled from cases involving travel from Tanzania, two of which are basal to the three sampled in Angola (Figure 1A, Supplemental Table S1). This new VOI, temporarily designated A.VOI.V2, has 31 amino acid substitutions (11 in spike) and three deletions (all in spike) (Figure 1C & 1D). The spike mutations include three substitutions in the receptor-binding domain (R346K, T478R and E484K); five substitutions and three deletions in the N-terminal domain, some of which are within the antigenic supersite (Y144{Delta}, R246M, SYL247-249{Delta} and W258L)4; and two substitutions adjacent to the S1/S2 cleavage site (H655Y and P681H). Several of these mutations are present in other VOCs/VOIs and are evolving under positive selection.

Freely accessible ready to use global infrastructure for SARS-CoV-2 monitoring (preprint)

The COVID-19 pandemic is the first global health crisis to occur in the age of big genomic data. Although data generation capacity is well established and sufficiently standardized, analytical capacity is not. To establish analytical capacity it is necessary to pull together global computational resources and deliver the best open source tools and analysis workflows within a ready to use, universally accessible resource. Such a resource should not be controlled by a single research group, institution, or country. Instead it should be maintained by a community of users and developers who ensure that the system remains operational and populated with current tools. A community is also essential for facilitating the types of discourse needed to establish best analytical practices. Bringing together public computational research infrastructure from the USA, Europe, and Australia, we developed a distributed data analysis platform that accomplishes these goals. It is immediately accessible to anyone in the world and is designed for the analysis of rapidly growing collections of deep sequencing datasets. We demonstrate its utility by detecting allelic variants in high-quality existing SARS-CoV-2 sequencing datasets and by continuous reanalysis of COG-UK data. All workflows, data, and documentation is available at https://covid19.galaxyproject.org.

Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa (preprint)

Continued uncontrolled transmission of the severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) in many parts of the world is creating the conditions for significant virus evolution. Here, we describe a new SARS-CoV-2 lineage (501Y.V2) characterised by eight lineage-defining mutations in the spike protein, including three at important residues in the receptor-binding domain (K417N, E484K and N501Y) that may have functional significance. This lineage emerged in South Africa after the first epidemic wave in a severely affected metropolitan area, Nelson Mandela Bay, located on the coast of the Eastern Cape Province. This lineage spread rapidly, becoming within weeks the dominant lineage in the Eastern Cape and Western Cape Provinces. Whilst the full significance of the mutations is yet to be determined, the genomic data, showing the rapid displacement of other lineages, suggest that this lineage may be associated with increased transmissibility.

Major new lineages of SARS-CoV-2 emerge and spread in South Africa during lockdown. (preprint)

In March 2020, the first cases of COVID-19 were reported in South Africa. The epidemic spread very fast despite an early and extreme lockdown and infected over 600,000 people, by far the highest number of infections in an African country. To rapidly understand the spread of SARS-CoV-2 in South Africa, we formed the Network for Genomics Surveillance in South Africa (NGS-SA). Here, we analyze 1,365 high quality whole genomes and identify 16 new lineages of SARS-CoV-2. Most of these unique lineages have mutations that are found hardly anywhere else in the world. We also show that three lineages spread widely in South Africa and contributed to ~42% of all of the infections in the country. This included the first identified C lineage of SARS-CoV-2, C.1, which has 16 mutations as compared with the original Wuhan sequence. C.1 was the most geographically widespread lineage in South Africa, causing infections in multiple provinces and in all of the eleven districts in KwaZulu-Natal (KZN), the most sampled province. Interestingly, the first South-African specific lineage, B.1.106, which was identified in April 2020, became extinct after nosocomial outbreaks were controlled. Our findings show that genomic surveillance can be implemented on a large scale in Africa to identify and control the spread of SARS-CoV-2.